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NUMERICAL DIGITAL COMPUTER METHOD FOR DETERMINING 
THE TRANSIENT RESPONSES OF NONLINEAR AUTOMATIC SYSTEMS 
BASED ON CALCULATION OF THE CONVOLUTION INTEGRAL 

by 

A. V. Vul'fson 


Discussed is a technique for the digital computer calculation of transient 
processes for systems with one or more nonlinear characteristics, using an 
extension of the convolution technique developed by Carson for systems with one 
nonlinearity. The method does not require formulation of a system of first-order 
differential equations with subsequent programing of the right-hand sides for 
each problem. The output data are transfer functions of the linear part of the 
system. The procedure for programing the solution of a specific problem is 
simplified, reducing essentially to the mere input of numerical data. The 
nonlinearities may be given tabularly, and they may be discontinuous. 


Transient processes of automatic systems are usually determined by 
solving numerically differential equations with the aid of digital computers . 
However, the investigation of a system with a complex linear part containing 
differentiating links, and of a system with one or several nonlinearities, 
encounters difficulties in reducing differential equations to a normal Cauchy 
form, particularly when the nonlinearities have breaks or discontinuities. 
Similar difficulties are also encountered in determining the transient processes 
of electrical circuits containing one or several nonlinear components. In many 
cases it is possible to avoid these difficulties by using the numerical methods 
of solving integral equations. 

J. Carson was probably the first to use a convolution integral in 
recording an equation for the transient process of a system with a branched 
linear part and one nonlinearity [1] . A numerical-graphic method of calculating 
the transient processes of similar systems, based on the calculation of the 
convolution integral, was developed by N. I. Sokolov [2]. 

A method is offered in this article of using a digital computer for 
determining the transient process of a system containing one or several 
nonlinear characteristics, which is based on the idea of Carson. This method 
requires no composing of systems of differential equations of the first order 
to be solved for the derivatives, and no programing for each problem of the 
right parts of the differential equations. The initial data consist of the transfer 
functions of the linear part of a system; it simplifies the programing procedure 
for solving a specific problem and reduces it mostly to an introduction of the 
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numerical material. The nonlinearities may be specified in tabular form and 
may have discontinuities (relay characteristics) . 

Numerical Determination of Transient Process 
of a System with One Nonlinearity 



FIGURE 1. A SYSTEM WITH FIGURE 2. REDUCED STRUCTURAL 
ARBITRARY STRUCTURE PATTERN OF A SYSTEM WITH 

OF LINEAR PART AND ONE NONLINEARITY 

WITH ONE NONLINEARITY 


The system shown in Figure 1 consists of a linear part with an arbitrary 
structure and one inertialess link with a static characteristic i = ip (u) (it is 
shown separately on the drawing) . The specified input action and the sought 
output response at the output of the system are, respectively, f(t) and x(t) ; 
u(t) and i(t) are, respectively, the input and output of the nonlinear link. If 
f(t) has a Laplace representation F(p) (it will be also assumed that this 
representation is a fractional-rational function) , the following system of 
equations can be written to determine x(t) [ 1] : 


X(P) = W (p)I(p) + 

X 

W xf (p)F(p), 

(a) 


U(P) = W(p)I(p) + 

W f (p)F(p), 

(b) 

(1) 

i(t) = !p [u ( t) ] . 


(c) 



At t < 0, it is assumed that x = u = i = 0. A similar writing corresponds 
to a system reduced to the form shown in Figure 2. The terms W (p) , 

X 

W xf (p) , W(P), and W^(p) are determined by the structure of the linear part of 

the system with the aid of the known methods of transforming the structures of 
linear systems. Let us assume that w (t) and w(t) are the originals of the 

X 

Laplace representations W (p) and W(p) 

X 
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f xf (t) = W xf (p)F(p), 

f f (t) = W f (p)F(p). 

Using the theorem of convolution, Equations (1) are represented in the 
form of a system containing integral equations [1, 2]: 

t 

x(t) = /w x (t - T)i(r)dr + f rf (t) , (a) 

0 

t 

u ( t) = /w(t - T)i(r)dT + f (t) , (b) (2) 

0 f 

i(t) = ip[u(t)}. (c) 


The first stage of the numerical calculation of the transient process 
consists of determining the transient characteristics w (t) , f f (t) , w(t) , and 

fj(t) . A numerical method is described [5] which can be used for composing 

a program which requires no changes for determining the transient character- 
istics from a fractional-rational Laplace representation of any order. This 
method is based on the use of interpolating expressions of Adams and provides 
a solution with an accuracy of the second order. 

It will be assumed that this method is used to determine the transient 
characteristics w (t) , f (t), w(t) , and f (t) . The latticed functions obtained 

X XI I 

as a result of the numerical calculation will be designated by w x [n] , f^n] , 

w[n], and f^[n]. Their quantization period is constant as is equal to the step 

At of the numerical determination of the transient characteristics. To find 
the numerical solution of the system ( 2) , the integrals entering ( 2a) and ( 2b) 
are replaced with finite sums of the quadrature formulas which are obtained as 
a result of the piecewise-linear interpolation of each of the discrete sequences: 
w [n] , w[n] , and i[n] . 

X 
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w [n - 1] w [n] 
x[n] = I + -V- I Ati[0] 


n-1 / w [n-m-1] w [n - m] w [n - m + 1] 

+ Z (— + -th-T + — 7 ) Ati[m] 

m = l ' 


1.5 


■) 


') 


/ w [0] w [1] 

- H— + -V) Atiini +f X f |n) ' 


(a) 


u[n] = 


w 


■ - 1 + Ati[0] 


+ 


n-1 

E 


m = l 


' w[n - m - 1 ] w[n - m] w[n - m 
L 6 + 1.5 + 6 


+ T) + V” 1 ' 

i[n] = $(u[n] ) . 
n = 1, 2, 3 . . . 




Ati[m] 

(b) 

(c) 


(3) 


Using the first differences of (3a) for (3b) , we reduce the system (3) 
to a form of: 


Ax[n] = R x [n] + S Ai[n], 


(a) 


Au[n] = R[n] + SAi[n], 

i[n] + Ai[n] = ^(u[n] + Au[n]) 


(b) (4) 

(c) 


n = 0, 1, 2 . . 


In which case. 


■( 

-( 


w [0] Aw [0] 

X X 


) At, 


wlOJ 


Aw [ Q 1 \ 


/ 


At: 


(5) 

(6) 
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/ Aw [0]\ 

R x t0] = l W x [0] + 2 j AU[0] + M xf t0] ’ 


( Aw [n - 1] Aw [n] \ 

—6 + Ati[ ° ] 


n - 1 / 
e = 1 ' 


(V) 


Aw [e — 1 ] Aw [e ] Aw [e + 1 ] \ 

^ — + -rr~ * —* — ) Au[n - c j 


/ 5 \ 

+ |w x [0] + -Aw x [0] + ) Ati[n] + Af rf [n]. 


n = 1, 2, 3 . . . 


R[0] = ^w[ 0] + A ^ [ ° ] ) Ati[0] + Af f [0], 

R[n] = (AwJ^] + Ati[0] 

■p"^) Ati[n - €] 

+ (w[0] + |aw[ 0] + Ati[n] + Af f [n]. 

n = 1, 2, 3 . . . 


^ / Aw [ € - 1 ] 


e = 1 


Aw[ € ] 
1.5 


Aw[ 


( 8 ) 


Note that for R [n] and R[n] the special expressions for n = 0 are 

X 

obtained only when the signal at the input of the nonlinearity or at the output of 
the system changes jumpwise when a perturbation is applied at the input of the 
system. In such a case, it is first necessary to determine: 


u[ 0] 

il 

© 

(a) 


i [ 0] 

= ^(u[ 0] ) , 

(b) 

(9) 

x[0] 

O 

II 

(c) 
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Therefore, the second stage of the numerical calculation of the transient 
process consists of determining for the nth step (n = 0, 1, . . . ) the value of 
R [n] and R[n] and solving the system of equations (4) for the unknown Ax[n], 

Au[n] , Ai[n], followed by the determination of: 


x[n + 1] = x[n] + Ax[n], 

(a) 


u[n + 1] = u[n] + Au[n], 

(b) 

(10) 

i[n + 1] = i[n] + Ai[n] . 

(c) 



The latticed function with a quantization period of At - x[n] is the one 
which represents the sought response at the output; if necessary, it can be 
interpolated, in which case it is expedient to use a parabolic interpolation when 
taking into account the order of accuracy of the expressions for the numerical 
determination of the transient characteristics and the determination of the 
convolution. 

It was stated above that the action at the input f(t) should have a 
fractional-rational Laplace representation. 

In those cases where this is not observed (for example, when f(t) is 
specified in tabular form) , f x ^ [ n. J and f^[n] can be determined by using the 

expressions for numerical determination of the convolution integral, as in 
Eq. (3) . In any other respects, the pattern of the solution remains unchanged. 

Conditions under Which the Method Can Be Used 

1. Each individual transfer function W (p), W(p) , W (p), and W (p) 

X XI I 

should correspond to a stable system. In other words, the linear system 
formed by opening a nonlinear link should be stable. An ultimate case of still 
being able to use the method is the presence of an integrating factor 1/p in the 
transfer functions. In such a case, a stable system should correspond to the 
transfer function remaining after separation of the integrating factor. 

2. The transfer functions w (t) , w(t) , f (t) , and f (t) should be finite 

X XI I 

when t = 0. Note that, as a rule, this requirement is satisfied by automatic 
control systems and, by no means always, by electrical circuits with passive 
nonlinear components. 
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The subsequent determination of the numerical solution is similar to 
the one used for a case with one nonlinearity. The only thing to be noted is that, 
instead of calculating four transient characteristics in case of one nonlinearity, 
it is necessary to calculate (d + l) 2 of such characteristics. A system of 
algebraic equations similar to (4) , in case of d nonlinearities, will be written 
as: 


Ax[n] = R x [n] + S^AiJn] + S x2 Ai 2 [n] + ... 
Au! [n] = Ri [n] + SjjAi^n] + S 12 Ai 2 [n] + ... 

Au 2 [n] = R 2 [n] + S 21 Ai 1 [n] + S 22 Ai 2 [n] + ... 

Au d [n]= R d [n] + S^AiJn] + S d2 Ai 2 [n] + ... 


+ 


+ 


+ 


+ 


xd 


Id 


2d 


dd 


Ai^n], 

^i d [n], 

Ai d [n], 

Ai d fn], 


( 12 ) 


ii [n] + Aij[n] = ip t (u t [n] + Au t [n] ) , 
i 2 [n] + Ai 2 [n] = ip 2 (u 2 [ n] + Au 2 [n]), 


i d [n] + Ai d [n] = ^(u d [n] + Au d [n] ). 
n = 0, 1, 2 . . . 

Here, remaining unknown, will be: 

Ax[n], Au t [n], Au 2 [n] Au^n], Ai t [n], Ai 2 [n], .... Ai d [n]. 

It is obvious that the time-and-labor spent for the solution increases propor- 
tionally to the square of the number of nonlinearities. 

The Experiment 

This method was used for programing the "Ural-2" digital computer. 
The program served to determine the transient process of a system containing 
one nonlinearity and a linear part having an arbitrary structure. The prepara- 
tion for the solution consists of composing Equations (1) for the system. This 
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Estimating the Accuracy of the Solution 


The most practically convenient method is the well-known method of the 
reduced step of calculation. Since the step remains constant during the solution 
of a problem, the problem should be solved twice and, during the second time, 
the step should be reduced to one-half, for example. 


In estimating the accumulated error, it should be assumed that the error 
is decreasing proportionally to the square of the step. It is expedient to estimate 
the accuracy only when the closed system is stable. 


Numerical Determination of Transient Processes of Systems 
Containing Several Nonlinear Links 



Under consideration is a system 
with d nonlinearities (their static 
characteristics are: i t = ^(u^ , 
i 2 = h( u 2 ) . • • • i d = ^ d ( u d > and whose 

linear part has an arbitrary structure 
( Figure 3) . As it was done for a case 
with one nonlinearity (1) , let us write 
for the transient process of the system 
of equations: 


FIGURE 3. A SYSTEM 
WITH A LINEAR 
PART HAVING AN 
ARBITRARY 
STRUCTURE AND 
CONTAININ' u 
NON LINE AF TIES 


X(p) = W^pJMp) + W x2 (p)I 2 (p) + ... 

+ w xd (p)! d (p) + W xf (p)F(p), 

Uj(p)= WnfplMp) + W 12 (p)I 2 (p) + ... 

+ W ld (p) I d (p) + W lf (p)F(p), 

U 2 (p)= W 21 (p)I 1 (p) + W 22 (p)I 2 (p) + ... 

+ W 2d (p) I d (p) + W 2f (p)F(p), 


U d (p) = W dl <fWP> + W d2 (p)I 2 (p) + ... + w dd (P)yP> + W df (p)F(p), 

ii(t) = ^ ! [uj (t) ] , 

Mt) = 4>2[u 2 ( t)], 





can be accomplished by using any of the methods of composing the transfer 
functions of linear systems (the method of directional graphs, for example) . 
Since following the reduction to the form of system (1) the algorithm for the 
solution is the same for any problem, no additional programing is required and 
the preparation for the counting is reduced to an introduction of numerical 
information: the numerator and denominator coefficients of the fractional- 
rational representations of W x ( p) , W xf (p)F(p), W(p), W f (p)F(p), theirorder, 

their step, and number of steps. Additional programing of i = ^ (u) is required 
only when a nonlinearity is specified analytically. When specified in tabular 
form, only the following numbers are introduced for the characteristics of the 
nonlinearities: the constant step u and the ordinates i. The program performs 
the sampling of the intermediate values by a parabolic interpolation. The 
system of algebraic equations ( 4) is solved by the method of iteration. 

As a result of the performance by the program, we obtain: 
x[n], u[n] , i[n], n = 0, 1, 2, . . . 

The program was used to determine the transient processes of several 
automatic systems (including relay systems) and of electrical circuits. In 
cases of periodic steady-state conditions, the fluctuations (forced or natural) 
were obtained merely by continuing the counting of the transient process. 

It was established that in case of linear parts having complex structures, 
only the composing of the system (1) becomes more time-consuming. The total 
time spent by the machine in solving the problem practically does not increase 
when the orders of the transfer functions in (1) increase; it is determined 
only by the number of steps. 

Let us consider one of the solved examples. Shown in Figure 4 is a 
block-diagram of a primary astatic control of the absolute angle of a turbine 
unit connected with a high-power system, during the correction through the 
turbine [4] . The reduced transfer function of the regulator is: 


W (p) 
P 


0. 00915p 5 + 0. 2305p 4 + 2. 4319p 3 + 14. 21p 2 + 30. 575p + 15. 03 
(0. 00545p 4 + 0. 1333p s + 1.581p 2 + 10.19p+ 22. 6) p 


The transfer function of the turbine unit is 


W Q (p) 


1 

0. 019p 2 + 0. Olp ‘ 
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LEGEND: 



FIGURE 4. BLOCK-DIAGRAM 
OF A SYSTEM OF 
PRIMARY ASTATIC 
CONTROL OF ABSOLUTE 
ANGLE GF A TURBINE 
UNIT CONNECTED WITH A 
HIGH-POWER SYSTEM 


<p = the deviation of the absolute angle 
jii = the increment of the power -exchange 
between the system and the unit 
H = the increment in power of the working sub- 
stance at the input of the turbine 

M = the load increment. 

11 


The system is nonlinear, because the 
relationship for the increment in power- 
exchange between the unit and the system and 
the deviation of the absolute angle is 
specified by the expression 

H 0 = 0.3 sin < p . 

Under investigation is the transient 
process in the system when a nominal load 
iUjj = 1 is added. 



FIGURE 5. CURVE OF THE 
TRANSIENT PROCESS 
cp( t) WITH LOAD ADDED 


The system corresponding to 
Equations (1) will be written for 
this case (in Laplace representa- 
tion) , as follows: 


W W 

= " 1 + W W P o (P ^ 
o p 


W W 

° P • I ( o) 

1 + W W p ’ ^ ; 
o p 


W 


<P(P) = ~ 


1 + W W 'o 


At (P) 


o p 


w 


O 1 

W W p’ 

op 

(b) 

0. 3 sin cpt. 

(c) 


( 13 ) 
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The selected step of solution was \t= 0. 04 seconds; the number of steps 
was 220. The time spent in introducing the information into the machine, after 
it was reduced to the form of system (13) , ranged from 25 to 30 minutes. The 
problem was solved by the machine in 16 minutes (including the control) . The 
curve of the transient process cp(t) is shown in Figure 5. The accuracy of the 
solution (in percent of maximum amplitude) was 0.65 percent. 


11 


preceding page blank not filmed. 


LITERATURE CITED 


1. J. R. Carson, THEORY AND CALCULATION OF VARIABLE ELECTRICAL 
SYSTEMS, Physical Review , No. 2, 1921. 

2. N. I. Sokolov, ANALYTICAL METHOD OF APPROXIMATE CALCULATION 
OF TRANSIENT PROCESSES OF CERTAIN NONLINEAR AUTOMATIC 
CONTROL SYSTEMS, Trudy MAI (Trans, of Moscow Aviation Institute) , 

No. 12, 1959. 

3. I. S. Berezin and N. P. Zhidkov, ME TODY VYCHISLENIY (Methods of 
Computing), 1-2 , Fizmatgiz, 1960. 

4. Ye. I. Yurevich, DEVELOPMENT OF AUTOMATIC CONTROL SYSTEM 
OF SUPERPOWER COMBINED SYSTEMS BY CONTROLLING THE ANGLE, 
Dissertation for the Scientific Degree of Doctor of Technical Sciences, 
Leningrad, 1963. 

5. A. V. Vul'fson, NUMERICAL METHOD OF DETERMINING THE TRANSIENT 
PROCESSES OF LINEAR AUTOMATIC SYSTEMS AND THEIR CALCULA- 
TION BY DIGITAL COMPUTERS. 


13 



DISTRIBUTION 


No. of Copies 

EXTERNAL 


Air University Library 1 

ATTN: AUL3T 

Maxwell Air Force Base, Alabama 36112 

U. S. Aimy Electronics Proving Ground 1 

ATTN: Technical Library 

Fort Huachuca, Arizona 85613 

U. S. Naval Ordnance Test Station 1 

ATTN: Technical Library, Code 753 

China Lake, California 93555 

U. S. Naval Ordnance Laboratory 1 

ATTN : Library 

Corona, California 91720 

Lawrence Radiation Laboratory 1 


ATTN: Technical Information Division 

P. 0. Box 808 

Livermore, California 94550 


Sandia Corporation 1 

ATTN : Technical Library 

P. 0. Box 969 

Livermore, California 94551 

U. S. Naval Postgraduate School 1 

ATTN: Library 

Monterey, California 93940 

Electronic Warfare Laboratory, USAECOM 1 

Post Office Box 205 

Mountain View, California 94042 

Jet Propulsion Laboratory 2 

ATTN: Library (TDS) 

4800 Oak Grove Drive 
Pasadena, California 91103 

U. S. Naval Missile Center 1 

ATTN: Technical Library, Code N3022 

Point Mugu, California 93041 

U. S, Army Air Defense Command 1 

ATTN: ADSX 

Ent Air Force Base, Colorado 80912 

Central Intelligence Agency 4 

ATTN: OCR/ DD- Standard Distribution 

Washington, D. C. 20505 

Harry Diamond Laboratories 1 

ATTN : Library 

Washington, D. C. 20438 

Scientific 6 Tech. Information Div., NASA 1 
ATTN: ATS 

Washington, D. C. 20546 


14 


No. of Copies 


U. S. Atomic Energy Coirenission 1 

ATTN: Reports Library, Room G-017 

Washington, D. C. 20545 

U. S. Naval Research Laboratory 1 

ATTN: Code 2027 

Washington, D. C. 20390 

Weapons Systems Evaluation Group 1 

Washington, D. C. 20305 

John F. Kennedy Space Center, NASA 2 

ATTN: KSC Library, Documents Section 

Kennedy Space Center, Florida 32899 

APGC (PGBPS-12) 1 

Eglin Air Force Base, Florida 32542 

U. S. Army CDC Infantry Agency 1 

Fort Benning, Georgia 31905 

Argonne National Laboratory 1 

ATTN: Report Section 

9700 South Cass Avenue 
Argonne, Illinois 60440 

U. S. Aimy Weapons Command 1 

ATTN: AMSWE-RDR 

Rock Island, Illinois 61201 

Rock Island Arsenal 1 

ATTN: SWERI'RDI 

Rock Island, Illinois 61201 

U. S. Army Cmd. 5 General Staff College 1 

ATTN: Acquisitions, Library Division 

Fort Leavenworth, Kansas 66027 

Combined Arms Group, USACDC 1 

ATTN: Op. Res., P and P Div. 

Fort Leavenworth, Kansas 66027 

U. S. Army CDC Armor Agency 1 

Fort Knox, Kentucky 40121 

Michoud Assembly Facility, NASA 1 

ATTN: Library, I-MICH-OSD 

P. 0. Box 29300 

New Orleans, Louisiana 70129 

Aberdeen Proving Ground 1 

ATTN: Technical Library, Bldg. 313 


Aberdeen Proving Ground, Maryland 21005 

NASA Sci. 6 Tech. Information Facility 5 

ATTN: Acquisitions Branch (S-AK/DL) 

P. O. Box 33 

College Park, Maryland 20740 

U. S. Army Edgewood Arsenal 1 

ATTN: Librarian, Tech. Info. Div. 

Edgewood Arsenal, Maryland 21010 



No. of Copies 


National Security Agency 1 

ATTN: C3/TDL 

Fort Meade, Maryland 20755 

Goddard Space Flight Center, NASA 1 

ATTN: Library, documents Section 

Greenbelt, Maryland 20771 

U. S. Naval Propellant Plant 1 

ATTN : Technical Library 

Indian Head, Maryland 20640 

U. S. Naval Ordnance Laboratory 1 

ATTN: Librarian, Eva Liberman 

Silver Spring, Maryland 20910 

Air Force Cambridge Research Labs. 1 


L. G. Hanscom Field 
ATTN: CRMXLR/Stop 29 

Bedford, Massachusetts 01730 


Springfield Armory 1 

ATTN: SWESP-RE 

Springfield, Massachusetts 01101 

U. S. Army Materials Research Agency 1 

ATTN: AMXMR-ATL 

Watertown, Massachusetts 02172 

Strategic Air Command (OAI) 1 

Offutt Air Force Base, Nebraska 68113 

Picatinny Arsenal, USAMUCOM 1 

ATTN: SMUPA-VA6 

Dover, New Jersey 07801 

U. S. Army Electronics Command 1 

ATTN: AMSEL-CB 

Fort Monmouth, New Jersey 07703 

Sandia Corporation 1 

ATTN : Technical Library 

P. 0. Box 5800 

Albuquerque, New Mexico 87115 

ORA(RRRT) 1 

Holloman Air Force Base, New Mexico 88330 

Los Alamos Scientific Laboratory 1 


ATTN : Report Library 

P. 0. Box 1663 

Los Alamos, New Mexico 87544 

White Sands Missile Range 1 

ATTN: Technical Library 

White Sands, New Mexico 88002 

Rome Air Development Center (EMLAL-1) 1 

ATTN: Documents Library 

Griffiss Air Force Base, New York 13440 


No. of Copies 

Brookhaven National Laboratory 1 

Technical Information Division 
ATTN: Classified Documents Group 

Upton, Long Island, New York 11973 


Watervliet Arsenal 1 

ATTN: SWEWV-RD 

Watervliet, New York 12189 

U. S. Army Research Office (ARO-D) 1 

ATIN: CRD- AA- IP 

Box CM, Duke Station 
Durham, North Carolina 27706 

Lewis Research Center, NASA 1 

ATTN: Library 

21000 Brookpark Road 
Cleveland, Ohio 44135 

Systems Engineering Group (RTD) 1 

ATTN: SEPIR 

WTight-PatteTSon Air Force Base, Ohio 45433 

U. S. Army Artillery § Missile School 1 

ATTN: Guided Missile Department 

Fort Sill, Oklahoma 73503 

U. S. Army CDC Artillery Agency 1 

ATIN : Library 

Fort Sill, Oklahoma 73504 

U. S. Army War College 1 

ATIN: Library 

Carlisle Barracks, Pennsylvania 17013 

U. S. Naval Air Development Center 1 

ATTN : Technical Library 

Johnsville, Warminster, Pennsylvania 18974 

Frankford Arsenal 1 

ATTN: C- 2500 -Library 

Philadelphia, Pennsylvania 19137 

Div. of Technical Information Ext. , USAEC 1 

P. 0. Box 62 

Oak Ridge, Tennessee 37830 

Oak Ridge National Laboratory 1 

ATIN : Central Files 

P. 0. Box X 

Oak Ridge, Tennessee 37830 

Air Defense Agency, USACDC 1 

ATTN: Library 

Fort Bliss, Texas 79916 

U. S. Army Air Defense School 1 

ATIN: AKBAAS-DR-R 

Fort Bliss, Texas 79906 


15 


U. S. Army CDC Nuclear Group 
Fort Bliss, Texas 79916 


No. of Copies 
1 


INTERNAL 


No. of Copies 


Manned Spacecraft Center, NASA 1 

ATTN: Technical Library, Code BM6 

Houston, Texas 77058 

Defense Documentation Center 20 

Cameron Station 
Alexandria, Virginia 22314 

U. S. Aimy Research Office 1 

ATTN: STINFO Division 

3045 Columbia Pike 
Arlington, Virginia 22204 

U. S. Naval Weapons Laboratory 1 

ATTN : Technical Library 

Dahlgren, Virginia 22448 

U. S. Army Engineer Res. 6 Dev. Labs. 2 

ATTN: Scientific $ Technical Info. Br. 

Fort Belvoir, Virginia 22060 

Langley Research Center, NASA 1 

ATTN: Library, MS- 185 

Hampton, Virginia 23365 

Research Analysis Corporation 1 

ATTN: Libraty 

McLean, Virginia 22101 

U. S. Army Tank Automotive Center 1 

ATTN: SM0TA-RTS.1 

Warren, Michigan 48090 

Hughes Aircraft Company 1 

Electronic Properties Information Center 
Florence Ave. 5 Teale St. 

Culver City, California 90230 

Atomics International, Div. of NAA 1 

Liquid Metals Information Center 
P. 0. Box 309 

Canoga Park, California 91305 

Foreign Technology Division 1 

ATTN : Library 


Wright -Patterson Air Force Base, Ohio 45400 

Clearinghouse for Federal Scientific and 1 

Technical Information 
U. S. Department of Commerce 
Springfield, Virginia 22151 

Foreign Science 5 Technology Center, USAMC 3 
ATTN: Mr. Shapiro 

Washington, D. C. 20315 

National Aeronautics § Space Administration 2 
Code USS-T (Translation Section) 

Washington, D. C. 20546 


Headquarters 

U. S. Army Missile Conmand 
Redstone Arsenal, Alabama 35809 
ATTN: AMSMI-D 1 

AMSMI-XE, Mr. Lowers 1 

AMSMI-XS, Dr. Carter 1 

ANBMI -Y 1 

AM5MI-R, Mr. McDaniel 1 

AMSMI-RAP 1 

AMSMI-RBLD 10 

USACDC-LnO 1 

AMSMI-RBT 8 

AMSMI-RB, Mr. Croxton 1 

National Aeronautics 8 Space Administration 
Marshall Space Flight Center 
Redstone Arsenal, Alabama 35809 
ATTN: MS-T, Mr. Wiggins 5 

R-COMP-T, Mr. Trauboth 1 


16 


UNCLASSIFIED 


Security Classification 


DOCUMENT CONTROL DATA - R & D 

(Security classification ot title, body ot abstract and indexing annotation must be entered when the overall report is classified) 

I* ORIGINATING activity ( Corporate author) 

Redstone Scientific Information Center 
Research and Development Directorate 
U. S. Army Missile Command 
Redstone Arsenal, Alabama 35809 

24. REPORT SECURITY CLASSIFICATION 

Unclassified 

2b. GROUP 

N/A 

3. REPORT TITLE 

NUMERICAL DIGITAL COMPUTER METHOD FOR DETERMINING THE TRANSIENT 
RESPONSES OF NONLINEAR AUTOMATIC SYSTEMS BASED ON CALCULATION OF THE 

Sve^iyaH^ssSi^^^Seljnvkh Zavedeniv, Elektromekhanika, No. 8, 841-848 (1965) 

4. descriptive notes (Type of report and inclusive dates) 

Translated from the Russian 

5- AUTHORfS) (First name, middle Initial , last name) 




A. V. Vul'fson 




«• REPORT DATE 

9 August 1967 

74. TOTAL NO. OF PAGES 

19 

7b. NO. OF REFS 

5 

34. CONTRACT OR GRANT NO- 

N/A 

b. PROJEC T NO. 

N/A 

M. ORIGINATOR'S 

RSIC-697 

REPORT NUMBER(S) 

C. 

<L 

OTHER REPORT NO(S) (Any other numbers that may be assigned 
this report) 

AD 

10. DISTRIBUTION STATEMENT 




Distribution of this document is unlimited. 




11* SUPPLEMENTARY NOTES 

12. SPONSORING MILI TARY ACTIVITY 1 

None 

Same as No. 1 



»3. abstract 


Discussed is a technique for the digital computer calculation of transient processes 
for systems with one or more nonlinear characteristics, using an extension of the convolution 
technique developed by Carson for systems with one nonlinearity. The method does not require 
formulation of a system of first-order differential equations with subsequent programing of 
the right-hand sides for each problem. The output data are transfer functions of the linear 
part of the system. The procedure for programing the solution of a specific problem is 
simplified, reducing essentially^ to the mere input of numerical data. The nonlinearities may 
be given tabular ly, and they may be discontinuous. 


DD 


POfltM 

« MOV «• 


1473 


RBPLACSS DO FORM 1473, 1 JAN B4. 
OBSOLETE FOR ARMY US E. 


UNCLASSIFIED 

^ Security CUssiftcstion 


17 








UNCLASSIFIED 



UNC LASSIFIED 


Stcurity Cla s • 1 flea tlon 


18 








